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SUMMARY 


Calculations made by the Laplace Transform method are presented showing the 
response to abrupt longitudinal control motion of the Shuttle and of four 
other airplanes, two of conventional and two of delta-wing configuration. Of 
these airplanes, three have been used in power-off landing experiments at the 
Dryden Flight Research Center. The effects of airplane type, center- of- gravity 
location, pitch damping, lift coefficient, and position of the cockpit on the 
response in height and pitching velocity are presented. 

Of the airplanes studied, only the Shuttle has been found to experience 
difficulties due to lack of precision in controlling the flight path or a 
tendency for longitudinal pilot-induced oscillations. The results show that 
the conventional airplanes experience relatively small initial reversal of 
response in altitude of the center of gravity compared to the delta-wing 
airplanes, and of these, the Shuttle has considerably larger reversed response. 
Factors contributing to this difference are discussed. The control difficul- 
ties of the Shuttle cannot necessarily be attributed to the response reversal, 
however, because of the presence of other adverse effects due to time delay and 
rate limiting in the control system. 


INTRODUCTION 


In the past, designers of unusually large or heavy airplanes have been concerned 
with the lag in response of these airplanes during the landing approach and 
flare. A number of studies of this problem have been made, such as those 
given in references 1-U. In reference 1, which considered transport airplanes 
of conventional configuration up to a gross weight of about 1.33 x 10° N 
(300 000 lb), the increased lag was predicted to be not of serious concern. 
Later reports, references 2-4, considered airplanes of larger gross weight 
and slender delta-wing designs with eleven control. Longitudinal control of 
the delta-wing configuration was predicted to- be a cause of concern because 
of the initially incorrect response in height due to lift produced by the 
elevens. In addition, the effect of the large distance of the pilot from the 
landing wheels in very large airplanes was thought likely to cause increased 
dispersion in touchdown point and sink rate in landing. The use of direct 
lift control was suggested as a means of alleviating some of these problems 
(ref. 5). 
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Since the publication of these reports, mmerous large transport airplanes, 
including a delta-wing configuration, have been placed in service without 
encountering any serious operational difficulties. In some cases, direct 
lift control was considered in the design but was found to be unnecessary 
in practice. 

The problem of the adverse initial response of a delta-wing configuration has 
received increased interest in the case of the Shuttle because it was thought 
to contribute to a pilot-induced oscillation problem encountered in the final 
landing of five landings made in the approach and landing tests. In discus- 
sion of this problem, data on the height response characteristics of different 
types of airplanes did not appear to be readily available. Perhaps one reason 
for this lack of data is that most flight measurements of transient responses 
are made with internal instrumentation which record linear accelerations and 
angular rates, but not changes in altitude. Changes in altitude are difficult 
to obtain directly with sufficient accuracy and sensitivity during rapid 
maneuvers. In addition, such changes in altitude are not readily apparent 
to the pilot except possibly during the landing flare where the fixed ground 
reference is available. 

The purpose of this report is to present calculated altitude variations and 
other longitudinal response quantities during abrupt control maneuvers for 
the Shuttle and for a number of types of airplanes. The effect of certain 
changes in stability characteristics, such as those due to center-of-gravity 
position and pitch damping, are also presented. 

All except one of the airplanes taken as examples have been used in power- 
off landing experiments at the Dryden Flight Research Center. These calcu- 
lations, therefore, serve as a basis of comparison with the Shuttle results. 
The one which has not been tested, a delta-wine bomber, was selected because 
it is a delta-wing airplane with weight similar to the Shuttle. 

The data presented may be of value for correlation with the known characteris- 
tics of airplanes similar to those analyzed. In addition, they may give a 
general appreciation of the range of altitude response characteristics 
encountered with different airplane types. A complete study of the effect of 
these characteristics on the longitudinal control problems during landing 
would require man-in-the-loop studies either by means of pilot modeling or by 
manned simulation investigations using ground-based or flight simulators. 
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SYMBOLS 


terms in formula ( 26 ) 
coefficients in transfer functions 
terms in formula ( 29 ) 

constants of integration 

■ I 

wing chord 

differential operator, — 

ds 

altitude, positive up 

d6- 

gain of pitch damper, d(qc/2V7 

ratio of radius of gyration in pitch to wing chord 
coefficient multiplying transfer function for Z/6 

e 

coefficient multiplying transfer function for D0/6' 

e 

radius of gyration in pitch 

ratio of tail length to wing chord 

ratio of distance between center of gravity and 
cockpit to chord length 

ratio of distance between center of gravity and 
instantaneous center of rotation to chord length 

pitching moment 

mass 

wing area 

nondimensional distance, chord lengths 
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time 

term in formula (25) 
airspeed 

vertical velocity, positive down 
vertical displacement or vertical force, 
angle of attack 

real roots of characteristic equation 
complex roots of characteristic equation 
characteristic determinant 
elevator angle 

downwash angle 

pitch angle, positive up 

relative density factor, 

’ pSc 

air density 

elevator effectiveness factor C /C 

Zi f Zi 

0 a, 
e t 

phase angle (formula 25) 
phase angle (formula ko) 


Subscripts 

w wing or wing-body 

t tail 

dot over a quantity indicates differentiation with respect to time 


positive down 
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stability derivatives are defined in accordance with the following examples. 
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ANALYSIS 


Derivation of Equations of Motion 

Because this analysis is concerned with rapid longitudinal maneuvers, the 
assumption of constant airspeed is made. The derivation of the short-period 
longitudinal equations has been given in mmierous previous reports and text- 
books, for example in reference 6. The derivation is repeated here for 
completeness, and to provide a clear presentation of the manner in which the 
results depend on the separate effects of the wing-body combination and the 
tail. This approach was also used in reference J. 

The longitudinal equations in dimensional form are as follows: 


on 




dZ 




K/ 






U V ^ ^ u. r ^Jd. 


(. 1 ) 


These equations are non-dimensionalized by making the following substitutions 


5 
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1 1 / 

c 

^ ^ ^ 

V dt 


Ky - 



fro 

e Sc 


( 2 ) 


In addition, the vertical force and moment are expressed in coefficient form 
in accordance with the usual conventions. 


C- 
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The equations then become 


~ °V % ^ % 

w '^e 

<;> KyP% = Cyr.^ + C„ + s Cm 

W ^ 


Note that in these equations the coefficients of tail moment and vertical 
force, ’ *^z ’ forth, are based on wing area and wing chord. 



The following approximations are made in calculating the angle of attack at 
the wing and the tail. 
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Furthermore, the lag operation e is expanded in a power series, and only 
the first-order term is retained. 


-LD 

e 



The expression for the angle of attack at the tail then becomes; 


.. 


oL^= OC - - ID) i- IDO 


( 5 ) 
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If the expressions (k) and (5) are substituted in (3) and coefficients of a , 
and so forth are collected, the equations may be written in the following 
form: 


® ( 6 ) 

^/j-Ky ' D & ~ ^ 


where 
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The stability derivatives of the basic airplane are considered to be composed 
of contributions from the wing-body combination and from the tail. This 
separation of effects makes it possible to estimate the derivatives for 
either a conventional or a tailless configuration. In the case of a tailless 
design, the effect of the tail contributions are omitted. Theories and 
empirical data such as those summarized in reference 8 may be used to 
estimate the values of the derivatives. With a conventional configuration, 
the tail contribution to the rotary derivatives may greatly exceed the 
contribution of the wing-body combination. In this case, the effect of the 
wing-body combination on the rotary derivatives may be estimated roughly or 
omitted. Formula (7) may be used, along with the expression 0-ii . ~ ~L^Z = 


to estimate the effects of the tail. The effect of the elevator may be 
calculated from the equations : 


Cr = r 


r 




<e 
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Transfer Functions for Longitudinal Response 

The transfer functions relating various non-dimensional response quantities 
to elevator motion may be foimd from formula (6). The transfer function for 
angle of attack, a, is: 




( 8 ) 


The transfer function for pitching velocity, D0, is 



The transfer function for non-dimensional normal acceleration, D(a - 6) is 




(10) 




The characteristic determinant, A , is 

o 


4 = 'S' 


+ i G”, - 7 Cn C. - KyQ. 

- -h ^ Cyy,^ 


( 11 ) 
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The non-dimensional vertical velocity and change in altitude may be obtained 
by integrating the non+dimensional normal acceleration DCa - 0) once and twice 
respectively. Thus 
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d fe J 


( 12 ) 
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The vertical motion of a point in the fuselage other than the center of 
gravity may be obtained by combining the effects of the vertical motion 
of the center of gravity and the pitching of the aircraft. For example, 
the non-dimensional normal acceleration at a point located at chord- 
lengths ahead of the center of gravity is given by the expression: 






'3 




(12a) 


In order to obtain dimensional values of the response quantities, the 
following relations may be used: 
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Effect of Pitch Damper 


Most modern airplanes have longitudinal stability augmentation in the form 

of pitch damping. A complete representation of such systems requires the 

consideration of instriiment and servo dynamics, washout filters, structural 

filters, etc. As a first approximation, however, the effect of a pitch 

damper may be considered as a change in the values of the stability derivatives 

C and C . The augmented value of these derivatives may be determined 
Z m 

q q 

from the formulas 


where 


r- 
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The use of pitch rate feedback affects only denominator terms in the transfer 

functions relating airplane response to control input. As shown by 

formula (9), the values of C and do not occur in the nimierator of 

m Z 

q q 

the transfer function for D0/6 . These terms do occur in the niunerator 
of the transfer function for a/6 , formula (8), and hence in D'{a - 6)/i5g’5 

formula 10. The combination of terms in which these quantities occur is: 
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If the modified derivatives given by (l4) are substituted in this expression, 
the resiilt is 






<£ 




X 


= c. 
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This result shows that the effect of rate feedback disappears from the 
numerator terms. Physically, this result can be understood by realizing that 
the relation between airplane response and elevator motion is the same 
whether the elevator is moved by the pilot or by a feedback signal. All 
response quantities, such as pitch rate and normal acceleration, therefore, 
have exactly the same relative values whether or not feedback is used. 


Calculation of Response to a Longitudinal Control Input 

The response to an arbitrary longitudinal control input may be obtained by 
numerical integration of the .equations of motion, or, in the case of simple 
inputs, by the use of Laplace transforms. In the present study, the response 
following an abrupt control motion is of interest. For this purpose, response 
to impulsive, step, or ramp inputs are examined. These responses may be 
obtained by the Laplace transform method. An application of the response to 
a ranq) input is the initial motion which results if the elevator control 
encoTonters rate limiting. The effect of such rate limiting, however, is to 
nullify the effect of any pitching velocity feedback. The stability deriva- 
tives of the basic airplane should, therefore, be used under conditions of 
rate limiting. 

Laplace transform tables usually present the transfer functions in a form with 
the coefficients of the highest power of the operator factored out . The 
transfer function of interest for the response of vertical position, Z, to 
control command takes the form 
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D^(D^fb,D t b^) 


(16) 
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Factoring out the coefficients, of the higher powers of D in equations (lO) 
and (ll) gives for the constants in equation (l6): 
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The denominator of equation (l6) may have real or complex roots. In the case 
of real roots, equation (l6) may be written 


^ =K 






pYo^ocXofr) 


( 22 ) 


The inverse transform of this expression, which may be found, for example, in 
reference 9 (page 3^+7, entry 2.2ii9) is 
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( 23 ) 


The expression gives the time response of Z to a unit impulsive longitudinal 
control input, because the inverse transform of an impulse is iinity. 

In case the denominator of equation (16) has complex roots, equation (16) may 
be written 


z 


p"f-a,Dra^ 


(2lt) 
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The inverse transform of this expression (ref. 10, pg. 237, entry 02.201) is 


Z 





i- fit f 



(25) 


where 


^ ^ P 

U = X~' (3 

^ = 2 tan 

' «■ 


ta r> 


(i('ZX~a, ) 

a-a^oc 


The response to a step input is obtained by multiplying the transfer function 
by the Laplace transform of a step, l/D, and taking the inverse transform 
of the result. Similarly, the response to a ramp is obtained by multiplying 
the transfer function by l/D^ and taking the inverse transform of the result. 
The time response of a system to a step input is, therefore, the integral of 
the response to an impulse, and the response to a ramp input is the second 
integral of the response to an impulse. The expressions for these transforms 
are of too high order to be found in most tables of Laplace transforms. As 
a result, the inverse transforms must be obtained from basic principles. In 
these cases, the simplest procedure for obtaining these time responses 
appears to be to take the first and second integrals, respectively, of the 
expressions for the time response to an impulse. 

For example, in the case of real roots, the expression for the response of 
the vertical position to an impulsive longitudinal input, equation (23), has 
the f orm 

/ -oct - /r , 

(26) 
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The response to a 
(26), is 


step longitudinal input , obtained by integrating equation 





r 


As , 



(27) 


The constant of integration, C , 
when t = 0. Hence, ^ 



A/ A:: 

'u T 


is obtained from the condition that Z = 0 


In the same manner, equation (27) may be integrated to obtain the response of 
vertical position to a ramp command. The result is: 



-ocir 

e -h -t,e f- 
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Where 
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In the case of complex roots, the equation for the response to an impulse, 
equation 25, has the form 


-ixr 


Z siv(et-'f-) 


(29) 
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The response to a step longitudinal input is : 
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The response to a ramp is : 

-d: 
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(31) 


The expressions for C and C^ are not written out because their values 
may be more readily obtained numerically in an actual example. 


The equations for the response in pitching velocity are as follows. As shown 
by equation (9)» the transfer function has the form 


— - fc / 

4 h,P th^ 


(32) 
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where 
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In the case that the denominator has real roots, equation (32) may he written 


= kA ■) 

‘S'e c^r^cvy^?/ r)J 


(35) 


The inverse transform of this expression, which may he found, for example, in 
reference 9 (ps- 338 j entry 1.107) is: 


I>6> = JC 


~yt 
y) <? ' 


l>0 


H' 




(36) 


This expression gives the resppnse of 1^(9 to an impulsive control input. 
Again, the response to a step and to a ramp control input may he obtained 
hy integrating this expression. The response to a step is 
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The response to a ramp is 


- KJ 



^ 7*' O 



(38) 


If pitching acceleration is wanted, it may be obtained by differentiating 
these equations. 

In the case the denominator has complex roots equation (29) may be written 


££ _ , r ^ % 

4 " 


(39) 


The inverse transform of this expression is (ref. 9, pg. 3^+2, entry 1.303) 


VB = K 


oe\ L ^ ^ ^ ^ > 
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where 


f := tan 




r= 




Again, the response to a step or ramp motion of the control may be obtained 
by integrating these expressions. The response to a step is: 




[jl%- ^ J 


-pcos{ptf^J f 


(i+1) 
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The response to a ramp is : 




•f 2 oC(Sco^{ Gt 


Examples 

Sketches of the airplanes used in the analysis are given in figure 1. The 
dimensions, weights, inertias, and longitudinal stability derivatives are 
given in table I. These derivatives were obtained partly from existing 
reports and partly from estimation by methods given in reference 8 or by use 
of equations (2) presented previously. All data for airplane 2 and part of 
the data for airplane 1+ were obtained from reference 11. In any casej the 
data should not be considered as design values approved by the manufacturer, 
but merely as rough estimates adequate for the type of comparisons made in 
this report. 

In all cases, the stability derivatives are considered to be constant, 
independent of angle of attack, lift coefficient, or dynamic pressure. In 
addition, ground effects on the derivatives are neglected. As a result of 
these assumptions, the flight path through the atmosphere following a 
disturbance is independent of airspeed, dynamic pressure, or ground proximity. 
Changes in time histories of plotted quantities at different values of 
lift coefficient are simply the result of the . change in time scale as a 
fimction of airspeed. 
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i in Calculations 


RESULTS AND DISCUSSION 


Simplifications Used in Analysis 

In presentation of the calculated data, the simplifying assumptions used 
should be kept in mind. In an actual airplane, true impulsive, step or 
ramp inputs cannot be obtained because of lags and rate limiting in the con- 
trol system and control actuators. For a more realistic representation of 
an actual control system a complete digital or analog simulation would be 
required. The data presented, however, can be used to survey some of the 
fundamental differences in response resulting from airplane configuration 
and from the other variables studied. 
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Presentation of Calculated Results 


Response of the Shuttle to impulsive, step, and ramp inputs of the elevator 
are shown in figure 2(a). Similar results with stability augmentation in the 
form of rate damping are shown in figure 2(b). In order to obtain a basis 
of comparison of the magnitude of the three types of inputs, the impulse 
was assumed to have, a magnitude of 1 radian second, or the area under the 
impulse is 1 radian times 1 second. The step input has 'a magnitude of . 

1 radian, and the ramp input a slope of 1 radian per second. As a result, . 
the area under the curve of elevator angle versus time for the step input 
equals that -under the impulse after 1 second, and the area under this curve 
for the ramp input equals that under the impiilse after 2 seconds. These 
inputs have a very large magnitude, greater than that attainable in practice. 
Because of linearity of the calculated responses, however, the curves may 
readily be scaled down. to correspond to inputs of any desired magnitude, 
such as, for example, 0.1 or 0.01 times those plotted. 

The curves of figure 2 and subsequent figures are plotted in dimensional form 
for a lift coefficient typical of that used in the landing approach prior 
to flare. In the case of the Shuttle, this lift coefficient was taken as 0.6. 
The results, for a pull-up, of the Shuttle show a marked do-wnward movement of 
the center of gravity prior to the subsequent response in the upward direction. 
This initial reversal of response occurs whether the input has the form of 
an impulse, step, or ramp. The main difference between these three types of 
inputs is a progressively increasing lag in obtaining the desired altitude 
responses in going from the impulse , to the ramp. 

The effect of stability augmentation in the form of rate damping is to reduce 
the magnitude of the initial reversed response and to reduce slightly the 
time to obtain the desired response. The rate feedback, also reduces the 
steady-state pitching velocity resulting from a step command. In practice, 
if rate feedback were used, the sensitivity, of the controller might be 
increased to retain the same, steady-state response. In this case, the main 
effect of stability augmentation on. the initial altitude response would be 
a slight decrease in the time required for the altitude to start to respond 
in the desired direction. 

The response of airplane 1, a heavy bomber airplane of conventional configura- 
tion, is sho-wn in figure 3. The most obvious difference between this case and 
the previous one is the great reduction in the initial reversed response of 
the motion of the center of gravity. 
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The altitude responses of the center of gravity for a step elevator (or 
elevon) motion for all of the airplanes studied are compared in figure U. 

No stability augmentation is assumed. The two airplanes of conventional 
configuration, airplanes 1 and 2, show very little reversed response at the 
start of the motion. The three delta-wing airplanes show reversed responses, 
but that of the Shuttle is much greater than those of the other two. The 
reason for this difference is discussed subsequently. 

A more detailed survey of the effects of stability augmentation on the Shuttle 
is given in figure 5. The pitching velocity feedback is varied from 0 to 
1.2hl degrees elevon per degree per second pitching velocity. This value 
is comparable to the feedback gain (GDQ in Shuttle nomenclature) used in 
the actual Shuttle, though the stabil^.ty augmentation system of the Shuttle 
contains additional features not simulated herein. As shown by this figure, 
the pitching velocity feedback causes the elevon angle to decrease following 
the initial step input. As a result, the response is more nearly similar 
to that for an impulsive input, and the time for the correct altitude response 
of the center of gravity is flight ly reduced. 

The effect of lift coefficient on the altitude response of t):ie Shuttle and 
of airplane 1 are shown in figure 6. As mentioned previously, the effects 
of lift coefficient are entirely the result of the change in time scale. 

In other words, the airplane flying more slowly takes longer to traverse 
the same path. 

The effects of shifting the center of gravity rearward from its assumed 
initial location are shown in figure T for the Shuttle and for airplane 2, 
both without stability augmentation. Larger increments of center of gravity 
movement are shown for airplane 2 because this airplane had more static 
stability to begin with. In either case, the effects of center of gravity 
position on the response during the first few seconds are practically 
negligible. This result is true even though the Shuttle was longitudinally 
unstable (center of gravity aft of the maneuver point) for the rearmost 
center of gravity location used. 

A comparison of the vertical motion at the cockpit location with that at 
the center of gravity for three airplanes is shown in figure 8. The effect 
of vertical location of the cockpit was neglected in making these calcula- 
tions. In the case of the Shuttle, the instantaneous center of rotation 
■following an impulsive elevon motion is somewhat ahead of the cockpit. The 
center of rotation is estimated to be 55-7 feet ahead of the center of 
gravity, whereas the cockpit is i+9-5 feet ahead of the center of gravity. 
Despite this difference, the vertical motion of the cockpit for the first 
0. T second of the motion is practically negligible, probably because the 
response to lift due to angle of attack of the wing is starting to take 
effect during this period. The difference in vertical motion between the 
cockpit and the center of gravity resulting from the pitch angle of the air- 
plane is sho'vm as a dashed line. 
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In the case of airplane 3, which is also a delta-wing configuration, the 
instantaneous center of rotation is aft of the cockpit and the cockpit 
responds in the correct direction from the start of the motion. The same 
result is shown for airplane 1, a conventional airplane configuration. In 
all three cases, the effect of the cockpit location ahead of the center of 
gravity is to reduce the lag in reaching a given change in altitude by about 
0.5 second. 


REASONS FOR DIFFERENCES IN TRANSIENT RESPONSE 


As shown by the precedine calculations , the airplanes with conventional 
configuration exhibit a relatively small initially reversed motion of the 
center of gravity compared to the delta-wing configuration, and of these, 
the Shuttle has considerably larger reversed response. As noted in reference 2, 
one factor influencing this behavior is the position of the instantaneous 
center of rotation of the aircraft following an abrupt elevator movement. 

If this center of rotation is far ahead of the center of gravity, initial 
reversed motion of the center of gravity is increased. 


The position of the instantaneous center of rotation may be obtained from 
formula 12a, which gives the transfer' ^function for vertical acceleration at 
any point along the fuselage. The initial acceleration, as shown by the 
initial value theorem (ref. 9» pg. 26 t) may be obtained by considering only 
the terms with the highest powers of the operator D in the transfer function. 
If these terms are substituted in formula 12a and the expression equated to 
zero, the distance in chord lengths 1 to the instantaneous center of 

X 

rotation m^ be obtained. The result is 




- £/( Ky Czg^ 


p 






C, 




In the last two terms of the denominator of this equation, the values C 

and C are frequently assumed equal to zero for delta-wing configurations. 

D-x. 

This assimiption is perhaps made for lack of better information. In the case 
of a conventional configuration, the combination of terms C + hC C 

^Doi“6e 

may be shown by substituting values from formula (j) to go to zero if only the 
tail contributions to the derivatives are considered. 
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In either case, then, these terms may be neglected. The expression for 1 
then reduces to 





The values of 1 for the airplanes studied are listed in table I. The 

distance in meters (ft) from the center of gravity to the cockpit and to the 
instantaneous center of rotation for these airplanes are -given in table II. 

If the values for the variables from table I are examined, the ratio /C 

Z o m r 

oe oe 

is seen to be quite similar for the three delta-wing configurations. The 
value of however, is much greater for the Shuttle. The greater radius of 

gyration in pitch is therefore an important reason for the increased reversed 
response of the Shuttle. This difference in inertia characteristics probably 
results from the unusual design layout of the Shuttle, which incorporates a 
large central cargo bay with most of the vehicle equipment in the nose and 
tail sections . The carsro bay is mostly empty in the loading condition 
considered. 

The instantaneous center of rotation determines the initial motion of the 
airplane, but the subsequent response cannot be expressed so simply. Never- 
theless, _ a large radius of gyration in pitch would be expected to reduce the 
rate of increase of angle of attack following a control deflection and, there- 
fore , delay the buildup of lift in the desired direction. 


DISCUSSION OF FLIGHT EXPERIENCE 


As mentioned previously, airplanes 1-3 have been used in power-off landing 
studies at the NASA Dryden Flight Research Center to study the power-off 
landing problems of airplanes such as the Shuttle and NASA research airplanes. 
In no case did these airplanes experience difficulties in obtaining the 
desired precision of touchdown conditions. Airnlane i* has not been tested 
in the power-off condition, but has been used in service with no comnlaints 
about its longitudinal characteristics in landing. In the case of the Shuttle, 
however, the control of altitude during the approach and flare has been 
described as lacking in predictability. A pilot-induced oscillation was 
experienced in the fifth landing of the approach and landing tests when an 
attempt was made to touch down at a specific point on the runway. Pilot- 
induced oscillations have also been experienced in flight simulation studies 
of a variable-stability airplane with characteristics adjusted to represent 
the Shuttle (ref. 12). 
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The Shuttle control system in the approach and landing tests had other 
response characteristics which might adversely affect the precision of 
control. These characteristics include a time delay between controller 
motion and elevon response which varied from 0.3 to 0.5 sec, and rate limiting 
of the elevens to 20 degrees per second. Thus, the control difficulties 
experienced cannot necessarily be attributed to the delay in altitude response 
caused by the lift due to elevon deflection. In unpublished simulation 
studies improvements in all three factors have been shown to be required to 
obtain desirable control characteristics. 

The presence of these additional adverse factors in the actual Shuttle control 
system prevents the use of the calculated response of the various airplanes 
(fig. to devise any simple criterion for allowable magnitude of the reversed 

response in altitude. Obviously, the path response of the Shuttle exhibits a 
more severe reversed response than any of the other airplanes studied. In 
the case of the Shuttle, analytical studies of the control characteristics 
using a pilot model, or manned simulation studies of the coniplete system using 
ground-based or flight simulators have proved promising in analyzing the 
cause of the difficulties and suggesting improvements. Such studies are 
beyond the scope of the present report. Some comments are offered, however, 
on the relative importance to the pilot of motion of the cockpit and motion 
of the center of gravity. 

The primary objective of a satisfactory landing is to touch down on the main 
landing gear with a desirably low sink rate. Inasmuch as the main landing 
gear is usually at about the same horizontal location as the center of gravity, 
the vertical motion of the center of gravity can be considered as representa- 
tive of that of the main landing gear. 

In some fixed-base simulation studies, a question has arisen as to whether the 
pilot should be given a simplified visual display representing the altitude 
of the cockpit or the altitude of the center of gravity. In cases in which 
the response has considerable lag, the reduction in lag associated with the 
forward location of the cockpit reduces the tendency for pilot-induced 
oscillations. Inasmuch as the pilot in the actual airplane is in a position 
to sense the vertical motion of the cockpit, the argument has been advanced 
that the simulation should represent the motion of the cockpit. In effect, 
this viewpoint holds that if the pilot can control the cockpit along a smooth 
flare path to the desired position at touchdown, the rest of the airplane will 
follow along behind and the wheels will touch down with a low sink rate. 

In actual flight, of course, the pilot can sense both normal acceleration and 
pitching acceleration, and the resulting velocities and positions in altitude 
and pitch. For this reason, the pilot can learn to sense the position of the 
wheels with respect to the ground whether he is in the nose of the airplane 
or in some other location. It is known that some fighter or racing type 
airplanes have been controlled successfully from a cockpit aft of the center 
of gravity. Thus, the position of the cockpit in the nose of an airplane 
such as the Shuttle may not alleviate all the control problems associated with 
lag and reversal of response at the center of gravity. If any control 
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corrections are required, the motion of the wheels will initially he 
incorrect, resulting in possibly increased vertical velocity at touchdown. 
The capability of the pilot to avoid pilot-induced oscillations may be 
alleviated by his ability to sense pitching acceleration, which is the 
same regardless of his location in the airplane. These argimients indicate 
that further research is required to determine the effect of pilot location 
on his ability to make satisfactory landings. In analytical and simulation 
studies of landing control problems, every effort should be made to include 
the effects of all cues which can be sensed by the pilot. The data of the 
present report may aid in evaluating the results of such studies and in 
comparing the results of tests on different types of airplanes. 


CONCLUDING REMARKS 


Results of calculations are presented showing the response to abrupt longi- 
tudinal control motion of the Shuttle and of four other airplanes, three of 
which have been used in power-off landing experiments at the Dryden Flight 
Research Center. The effects of airplane type, pitch damping, center-of 
gravity location, lift coefficient, and position of the cockpit are presented. 
These results are given mainly to provide a survey of response characteristics 
for different types of airplanes, including conventional and delta-wing 
configurations. 

Airplanes with the control surface aft of the center of gravity exhibit' an 
initially reversed response in altitude caused by the lift due to control 
deflections. The airplanes with conventional configuration exhibit a 
relatively small initially reversed motion of the center of gravity compared 
to the delta-wing configurations, and of these the Shuttle has considerably 
larger reversed response. Because of the presence of unusually large effeqtq 
of time lag and eleven rate limiting in the Shuttle control system, however , 
control difficulties experienced by the Shuttle cannot necessarily be 
attributed to the delay in altitude response caused by the lift due to 
eleven deflection. 
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Shuttle 


2 


3 


Airplane 


m, kg (slug) 

82309 (56U0) 

1.183 X 10^(8103) 

10885 ( 745 . 9 ) 

34602 (2371) 

68037 ( 4662 ) 

p, kg/m^(slug/ft^) 

1.139 (.00221) 

1.139 (. 00221 ) 

1.139 (. 00221 ) 

1.139 (. 00221 ) 

1.139 (. 00221 ) 

S, (ft^) 

249.9 (2690) 

371.6 ( 4000 ) 

18.22 (196.1) 

149.1 (1605) 

143.2 (1542) 

c, m (ft) 

12.06 (39.57) 

6.995 (22.95) 

2.904 (9.53) 

11.49 (37.70) 

11.03 (36.20) 

I , kgm^(slug ft^) 

V 

8.729 X 10^^ 

( 6.438 X 10 *^) 

1.302 X lo [ 

(9.60 X 10°) 

li 

8.813 X 10. 
(6.50 X 10 ) 

1.600 X 10^ 
(1.18 X 10°) 

1.417 X 10 ^, 
( 1.045 X 10 °) 

K 

y 

0.8539 

1.500 

0.9795 

0.5917 

0.4136 

y 

23.97 

. 39.94 

H80.6 

17.73 

37.79 

» per rad 
a 

- 2.70 

-3.75 

- 4.44 

-3.15 

- 2.84 

C , per rad 
m 
a 

-0.029 

-0.187 

-1.496 

1 

-0.089 

-0'.j142 

C , per rad 

Li 

q 

0.0 

- 3.66 

0.0 

0.0 

-1.80 

C , per rad 
m 

q 

-2.778 

-10.76 

-5.61 

- 1 . 4 o 

-5.91 

C , per rad 

^Dct 

0.0 

-1.28 

0.0 

0.0 ■ 

- 1 . 4 l 

C , per rad 

™Da 

o.,o 

-3.77 

-3.44 

0.0 

0.0 

C „ j per rad 
^ 6 e 

-0.956 

- 0.249 

-0.762 

-0.630 

-1.136 

C , per rad 

“ 6 e 

-0.495 

-0.732 

-1.50 

-0.351 

-0.437 

1 

X 

1 . 4 l 

0.765 

0.487 

0.628 

0.445 









TABLE II 


COMPARISON OF DISTANCE FROM CENTER OF GRAVITY TO COCKPIT WITH 
DISTANCE FROM CENTER OF GRAVITY TO INSTANTANEOUS CENTER OF ROTATION 


Airplane 

c.g. to Cockpit 
m ( f t ) 

c.g. to Center of Rotation 
m (ft) 

Shuttle 

15.1 (U9.5) 

17.0 (55.7) 

1 

IT.7 (58.1) 

(17.6) 

2 

6.3 ( 20 . 6 ) 

l.lt (li. 6 ) 

3 

15.2 (U9.8) 

7.2 (23.7) 

k 

11.6 (38.2) 

It. 9 (16.1) 
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Airplane 1 Airplane 2 



Airplane 3 Airplane 4 


Figure 1.- Sketches of airplanes used as examples, drawn to same scale 






Altitude of c.g. 



a. a 1 .0 2.0 3.0 H.0 


Time, sec 


Figure k.- Comparison of vertical motion of the center of gravity 
following a step elevator input of 1 radian for all 
airplanes studied. No stability augmentation. 
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Altitude of c. g. ± Altitude of eg. 





Shuttle, C|_=0.6 



“0 
0.02 
a 04 



Altitude 



04 


Figure -8.- 


Comparison of altitude response at co 
at c.g. following a step input 
The effect of pitch angle on the diff 
cochpit and c.g. is also shown. 
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